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Indirect Search for Dark Matter from the center of the Milky Way with the 

Fermi-Large Area Telescope 
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Vincenzo Vitale and Aldo Morselli, for the Fermi/LAT CoUaboration 
Istituto Nazionale di Fisica Nucleare, Sez. Roma Tor Vergata, Roma, Italy 

The Galactic Center region is expected to host the largest density of Dark Matter (DM) particles 
within the Milky Way. Then a relatively large gamma-ray signal would be expected from the 
possible DM particles annihilation (or decay). We are searching for the DM gamma-ray signal from 
the Galactic Center, which is also rich in bright discrete gamma-ray sources. Furthermore intense 
diffuse gamma-ray emission due to cosmic-ray interactions with interstellar gas and radiation is 
detected from the same direction. 

A preliminary analysis of the data, taken during the first 11 months of the Fermi satellite op- 
erations, is reported. The diffuse gamma-ray backgrounds and discrete sources, as we know them 
today, can account for the large majority of the detected gamma-ray emission from the Galactic 
Center. Nevertheless a residual emission is left, not accounted for by the above models. 

An improved model of the Galactic diffuse emission and a careful evaluation of new (possibly 
unresolved) sources (or source populations) will improve the sensitivity for a DM search. 



I. INDIRECT SEARCH FOR DARK MATTER 
WITH HIGH ENERGY 7 RAYS 

After the measurement of the orbital velocities of 
galaxies within clusters [ij, many other evidences 
point to the existence of a form of matter which is 
not coupled to the electromagnetic or to the strong 
interactions. Among these it is possible to list: the 
accurate measurements of the cosmic microwave back- 
ground (CMBR) 2]; the the galaxies rotation veloc- 
ities [3j],iJ]; the gravitational lensing 5|; the abun- 
dances of light elements ; and the large scale struc- 
tures [zl. With the measurements of the CMBR [§] 
the 23% of the Universe is estimated to be Cold Dark 
Matter (DM), while only 4% ordinary baryonic mat- 
ter. 

DM is coupled with gravitational interaction, and 
via this force to the ordinary matter. As a matter of 
fact N-body simulations predict that DM is arranged 
in large structures called halos ([§]). The DM density 
in a halo is large at the centre, and decreases with the 
distance. The density profile in the central part of the 
halo is not experimentally known and is studied by 
means of fits of N-body simulations or with analytic 
approaches [Toj . 

Non-gravitational DM couplings are studied with: 

• the direct search for DM scattering on ordinary 
matter; 

• the indirect study of DM annihilation via the 
secondary products, both charged and neutral 
(e"*" , p, d, V , ^ rays and lower frequency elec- 
tromagnetic radiation) . 

DM particles might produce gamma rays: 

• if the DM particles self-annihilate in pairs, pos- 
sibly following the scheme in Fig. 1. This yields 
a continuum gamma-ray emission, which is pro- 
duced by hadronization (or final state radiation) 



of the annihilation products and has a cut-off at 
the DM particle mass. The direct production 
of two gamma rays as annihilation products is 
suppressed in many models (10^'^ - 10^"* of the 
continuum); 

• if pseudo-stable DM particles decay in gammas. 
In this case the gamma-ray flux is proportional 
to the DM particle density. The DM decay con- 
stant is bound to be larger than 10^^ sec and is 
model dependent 

The gamma-ray flux from DM particle annihilation 
can be decomposed into a particle physics factor and 
an astrophysical one. The gamma-ray flux from anni- 
hilation can be written as: 
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where a is the annihilation cross-section, v is the rela- 



tive velocity of the DM particles 



is the number 



of gamma rays emitted per annihilation for each an- 
nihilation channel f, B/ is the branching ratio of the 
annihilation channel f and yiwiM p is the mass of the 
DM particle candidate. p(l) is the DM density pro- 
file, ie the density of DM as a function of the distance 
from the halo center. p(l)^ is integrated along the ob- 
server's line-of-sight {las). The gamma-ray flux from 
decay can be written as: 
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where F is the decay constant, is the number of 
gamma rays emitted for each decay and Mdm is the 
mass of the DM particle candidate, pdm is the DM 
density profile. 
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FIG. 1: gamma-ray production scheme for DM annihi- 
lation. Reproduced from 12]. If the DM particles self- 
annihilate and produce quarks, leptons and gauge bosons, 
then an indirect search is possible by searching for the 
secondary products. For weakly interacting massive parti- 
cles (WIMPs) of Maiorana type the heavy fermions pairs, 
such as bb,tt ,t'^t~ , are favoured as annihilation products. 
The annihilation into two gamma rays, for some models. 



is loop-suppressed, with a branching ratio of 10 



-10" 



Gamma-ray emission is expected after the heavy fermions 
hadronize. The gamma-ray emission energy spectrum is 
a continuum with curved shape (in log-log scale) and a 
sharp cut-off at the DM particle mass. 



The strong dependence of the emitted 7-ray flux on 
the DM density profile is evident in the annihilation 
case, smaller in the decay one. In fact in the decay 
scenario p, not its square, is integrated along the ob- 
server's line-of-sight (los, see Fig. 2). 

In order to search for a gamma-ray signal from DM 
annihilation (or decay), we observe the region from 
which we expect the largest signals or the best signal- 
to-noise ratio, where the noise is the gamma-ray back- 
ground. 

Our Galaxy could be embedded in a DM halo with 
the central density enhancement co-located with the 
Galactic Center. As the DM gamma-ray emission is 
expected to be a function of the DM density profile 
the Galactic Center is expected to be the brightest 
source of gamma rays from DM annihilation (see Fig. 
2 and 3) . It might be possible to detect the relatively 
large DM gamma emission ([l^, [13], [3), despite the 
bright gamma-ray emission from astrophysical pro- 
cesses, which is generated both by discrete sources and 
diffuse backgrounds, in the Galactic Center region 

A better ratio between the Dark and baryonic mat- 
ter, and then a weaker gamma-ray background, is fore- 
seen for the Dwarf Spheroidal Galaxies, or in DM sub- 
structures in the halo Il6l. 
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FIG. 2: Comparison of expected gamma-ray emission pro- 
file in the case of annihilating (proportional to the square 
of the density) and decaying (proportional to the den- 
sity) Dark. The angular profile of the gamma-ray signal 
is shown, as function of the angle 6 to the centre of the 
galaxy for a Navarro-Frenk-White (NFW) halo distribu- 
tion for decaying DM, solid (red) line, compared to the 
case of self-annihilating DM, dashed (blue) line. Both sig- 
nals have been normalised to their values at the Galactic 
poles, 9 = ±90°. Reproduced from [Tl|. 
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FIG. 3: Simulated gamma emission from an annihilating 
DM halo. The total surface brightness from all compo- 
nents together (the main halo -I- all resolved subhalos -|- 
all unresolved subhalos down to the free streaming limit) 
is shown. The surface brightness is in color code. Repro- 
duced from 



II. FERMI LARGE AREA TELESCOPE 

The Fermi Gamma-ray Space Telescope is a space 
observatory for the study of gamma-ray emission from 
astrophysical sources. It is equipped with two power- 
ful instruments: 

1. the Large Area Telescope (LAT), a pair- 
conversion telescope, composed of a preci- 
sion silicon-tracker (18 double-sided layers) and 
a calorimeter (8.6 radiation lengths vertical 
depth) , each of which consisting of a 4 x 4 array 
of 16 modules and surrounded by a segmented 
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anti-coincidence detector (ACD). LAT is oper- 
ated as a gamma-ray imager in the energy range 
between 20 MeV to over 300 GeV; 

2. the gamma-ray Bm'st Monitor (GBM), a detec- 
tor covering the 8 keV-40 MeV energy range, 
devoted to the study of the gamma-ray Bursts. 

Detailed descriptions of the Fermi observatory can be 
found in llTII and the LAT on-orbit cahbration is re- 
ported in [18l |. 

Studies of the gamma-ray sources in the energy 
range between hundreds of MeV and tens of GeV 
were performed by means of the SAS 2 and COS B 
sateUites in the years 1972 - 1982 and with the En- 
ergetic gamma-ray Experiment Telescope (EGRET) 
onboard of the Compton gamma-ray Observatory be- 
tween 1991 and 2000. EGRET detected 271 sources 
[l9l |. An half of the EGRET sources are unidentified, 
mainly because of the relatively large errors associ- 
ated with the source locations. The sources detected 
with EGRET are mainly members of two classes: Ac- 
tive Galactic Nuclei (AGNs) of FSRQ (Flat Spectrum 
Radio Quasars) and BL Lac types, with powerful rel- 
ativistic jets of plasma, and pulsars (spinning neu- 
tron stars, with powerful magnetic fields, capable of 
accelerating particles up to the high energy regime). 
EGRET was able to detect both discrete sources and 
diffuse gamma-ray emissions such as the extra-galactic 
gamma-ray background and the very intense Galactic 
component. 

The Large Area Telescope has an effective area five 
times larger, a much better angular resolution, and 
a sensitivity more than 10 times better than its pre- 
decessor EGRET. The Fermi observatory has several 
scientific objectives, which span many topics of as- 
trophysics and fundamental physics: (1) the detailed 
study of pulsars, AGNs, diffuse emissions and gamma- 
ray emission from nearby bodies; (2) the study of 
gamma-ray Bursts, up to GeV energies with the LAT 
and in the better studied keV- MeV range by means of 
the GBM; (3) the search for new classes of gamma-ray 
emitters; (4) the possible signals of new physics. The 
indirect search for DM particles and the investigation 
of their nature are major research topics for Fermi. 
Here we report an update of the indirect search for 
DM from the Galactic Center (GC). 



III. THE GALACTIC CENTER DURING THE 
FIRST YEAR OF FERMI 

Before the launch of the Fermi observatory, the sen- 
sitivity of the LAT instrument for the indirect search 
for DM was estimated in |l2|. In this publication 
many possible targets for the observation of gamma- 
ray emission from annihilating DM are considered and 
the possible limits on the velocity-averaged annihila- 
tion cross-section (<(tv> parameter) are given, as a 



function of the DM particle mass. For the Galac- 
tic Center, the considered gamma-ray background 
was the Galactic Diffuse emission, but not discrete 
bright gamma-ray sources. The results are obtained 
by means of analysis. The plot in Fig. 4 shows the 
constraints in <ctv> vs WIMP mass for an assumed 
pure bb annihilation channel. 

A likelihood analysis was instead performed in [23] , 
on simulations of both annihilating DM gamma-ray 
emission and gamma-ray background (both Galactic 
diffuse emission and the known discrete gamma-ray 
sources). With this model it was possible to conclude 
that the Fermi/LAT is sensitive to many benchmark 
realizations of a DM source, under the assumption 
that the simulation of the gamma-ray background is 
realistic. For example a DM particle with mass = 50 
GeV and a <av> = 3x lO^^cm'^s^^ would have pro- 
vided «280 events, within a ROI of 1.8° radius, in 30 
Ms. In the same time the simulated background would 
have generated ^20000 events, and a detection at 5a 
would have been possible, with a maximum likelihood 
analysis, which exploits both the spatial and the spec- 
tral information. Unfortunately the simulated back- 
ground did not reproduce the real complexity of the 
region. 
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FIG. 4: The Fermi gamma-ray Satellite sensitivity for DM 
indirect searches was investigated in For the GC the 
DM mass was considered between 10 and 1000 GeV, while 
the <(Tv> parameter between 0.5 and lOOxlO^^cm^ s~^ . 
The Galactic diffuse emission (both conventional and op- 
timized GALPROP models) was assumed as background, 
while the sources were not included as background, and a 
analysis was applied. In Fig 3 is reported the region of 
the <CTv> vs DM mass plane which can be sampled with 
Fermi/LAT in 5 years, with the observations of the Galac- 
tic Center. For this plot a pure bb annihilation channel 
and a NFW spatial distribution are assumed. 

The first release of results on the Galactic Cen- 
ter region obtained with the Fermi/LAT instrument. 
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was made with the Fermi Bright Source List pub- 
Ucation ([2l|, see Table 1). The Fermi/LAT in 
3 months produced a deep and well-resolved map 
of the 7-ray sky, thus also providing much bet- 
ter localization of the previous known gamma-ray 
sources. 205 sources were detected during the first 
3 months with statistical significance larger than 
10 a. The minimum flux for detection of this 
quality changed from the Galactic plane regions 
($(E>100 MeV)=0.4xlO-6cm-2s-i) to the Galac- 
tic poles ($(E>100 MeV)=0.05xlO-'5cm-2s-i), as 
a function of the intensity of the bright astrophysi- 
cal gamma-ra y b ackground (mainly the galactic com- 



ponent). In [21| the source OFGL J1746.0-29.00 is 
the source closest to the dynamical center of the of 
the Galaxy. OFGL J1746.0-2900 was reported to be 
marginally variable in 21], but this result was not 
confirmed with larger statistics. 

OFGL J1746.0-2900 is likely associated to 3EG 
J1746-2851, an unidentified EGRET source. A 
re-analysis of the EGRET data [2^ indicated 
for this source a steady integral flux <i>(>100 
MeV) = (118±73)xl0"^cm"2s"\ with a position dis- 
placed from the exact Galactic center toward positive 
Galactic longitudes. 

The TeV source in the Galactic Center (HESS 
J1745-290), a bright source above 100 GeV, is no- 
tably located near OFGL J 1746. 0-2900. The energy 
spectrum measured with H.E.S.S. [23| over the three 
years of data taking is compatible with both a power 
law spectrum with an exponential cut-off (photon in- 
dex of 2.10±0.04sfat±0.10sysf and a cut-off energy at 
15.7±3.4st(jt ±2.bsyst TeV) and a broken power law 
spectrum (spectral indices of 2.02±0.08stat± Q.l^syst 
and 2.63±0.14stat±0.10sj,st with a break energy at 
2.57±0.19stat± 0.44sj^st TeV). No significant flux vari- 
ation was found. The identification of the TeV 
gamma-ray source is still pending, although the DM 
annihilation hypothesis for the TeV emission is not 
favoured Furthermore the H.E.S.S. collabora- 

tion reported also weaker dif fuse emission from the 
central part of the Galactic Ridge [1^ . 



Source 
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Int. Flux 










(1<E<100 GeV) 


OFGL id 


n 


n 


n 


(lO-^cm-^s-i ) 


J1732.8-3135 


356.287 


0.920 


0.087 


3.89±0.33 


J1741.4-3046 


357.959 


-0.189 


0.197 


2.00±0.31 


J1746.0-2900 


359.988 


-0.111 


0.068 


7.92±0.47 



TABLE L Sources of the Fermi Bright Sources List, with - 
5° <1<5° and -5° <b<5°. The source identifier, the galac- 
tic longitude and latitude, the radius of the 95% source 
location confidence region and the integral flux between 1 
and 100 GeV are reported. The Bright Sources were de- 
tected with a statistical significance above 10a in 3months. 

During summer 2009 a first preliminary analysis of 



the Galactic Center data, for the indirect search for 
DM, was reported both at the TeVPA [265 and ICRC 
2009 [27] international conferences. For those results 
the likelihood analysis was used: 

• a Region-Of-Interest (ROI) of l°xl° was con- 
sidered; 

• the ROI was centred at the position RA = 
266.46°, Dec=-28.97°; 

• the Data taken during the first 8 months 
(8/2008-4/2009) have been used; 

• the events were selected to have energy between 
200 MeV and 40 GeV; 

• only events classified of " diffuse" class were se- 
lected for the analysis. 

For that 8 months analysis data were binned into a 
single spatial bin, then no spatial information was con- 
sidered. Likelihood analysis for the LAT data requires 
a model of the background (both diffuse emission and 
known sources). At that time a broken power law 
was assumed as the model of the gamma-ray emis- 
sion from this small ROI, as it was well-describing the 
observed data. The broken power-law spectrum was 
found to have integral flux between 100 MeV and 100 
GeV $(0.1 - lOOGel/) = (1.22±0.02) 10"^ cm'^s"!, 
first spectral index Fi= -1.38 ± 0.04, second spec- 
tral index ri= -2.60± 0.05, energy break at ^break 
— (1.6±0.1)GeV. Then a component representing the 
gamma-ray emission of annihilating DM was added to 
the likelihood analysis model. A benchmark DM real- 
ization with mass = 50 GeV and a pure hh annihilation 
channel spectrum was considered. 

With this model it was possible to put an upper 
limit of 2.43±0.02 lO""^ cm^^g-i (95% confidence 
level, above 100 MeV) to the integral gamma-ray flux 
from a DM annihilation source, which was assumed 
to have a Navarro-Frenk- White density profile, and 
the above reported benchmark mass of 50 GeV. This 
gamma-ray fiux upper limit implies a limit of 39.8 
X 10^^^ cm'^s"^ to the <tTv> parameter, if a pure 
bb annihilation channel is assumed. Nevertheless the 
analysis is to be considered very conservative, because 
no spatial information was used. 



IV. PRELIMINARY ANALYSIS OF FIRST 11 
MONTHS DATA 



The reported results were obtained with a binned 
likelihood analysis, performed by means of the tools 
developed by the Fermi/LAT collaboration (gtlike, 
from the Fermi analysis tools [1^). For this analy- 
sis: 
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FIG. 5: Spectra from the likelihood analysis of the 
Fermi/LAT data (number of counts vs reconstructed en- 
ergy). The likelihood analysis is the standard one used 
with the LAT data. The main analysis steps are: (1) to 
select data of high quality (selection cuts on events energy, 
zenith angle, reconstruction and classification quality); (2) 
to build a emission model of the region, based on the previ- 
ous knowledge and experimental evidence of new excesses 
with enough statistical significance; (3) to apply the like- 
lihood analysis to the data and the considered model. For 
each model component a fit of the free parameters and the 
computation the statistical significance is obtained. Here 
in the plot, from highest to the lowest: the black point are 
the observed data; the black line is the sum of all the com- 
ponents; the red line is the Galactic diffuse emission; the 
lower black line is the isotropic extragalactic; the other 
components are the sources detected. These results are 
preliminary. 
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FIG. 6: Residuals ( (exp.data - model) /model ) of the 
above likelihood analysis. For this residual spectrum the 
model and data were integrated over the 7°x7° ROI. The 
residual spectrum suggests that an unmodeled excess is 
present in the ~2-5 GeV range. The blue area shows the 
systematic errors on the effective area. These results are 
preliminary. 



• a ROI of 7° X 7° was considered. This ROI was 
used in order to minimize the background con- 
tribution and to avoid significant leakage of the 
gamma-ray signal under study; 

• the ROI was centred at the position RA = 
266.46°, Dec=-28.97°, ie the position of the 



brightest source; 

• the Data taken during the first 11 months 
(8/2008-7/2009) have been used; 

• the events were selected to have energy between 
400 MeV and 100 GeV; 

• only events classified of diffuse class and which 
converted in the front part of the tracker have 
been selected for the analysis. The selection 
in energy, event classification and conversion 
provided us with events with very well recon- 
structed incoming direction and data have been 
binned into a 100x100 bins map; 

• the IRF and the events classification are 
those relative to the Pass6V3 version of the 
Fermi/LAT analysis software. 

In order to perform the likelihood analysis for the 
LAT data, a model of the already known sources and 
the diffuse background should be built. The model in 
use for the presented analysis contains 11 sources in 
the Fermi 1 year catalog (to be published) which are 
located into the considered ROI, or located very close 
to the ROI boundaries and have a significant fraction 
of their flux leaking within the studied region. These 
sources have a point-like spatial model and a spectrum 
in the form of a power-law. The model also contain 
the diffuse gamma-ray background which is composed 
by two components: 

1. the Galactic Diffuse gamma-ray background. 
The observed Galactic Diffuse emission was 
modelled by means of the GALPROP code 
and [lOl, and the realization of the galactic 
emission named gll_iem_54_87Xexph7S.fit was 
used. During the likelihood maximization only 
the normalization of the GALPROP model is 
varied, not its components; 

2. the Isotropic Background. This component 
should account for both the Extragalactic 
gamma-ray emission and residual charged parti- 
cles. It is modelled as an isotropic emission with 
a template spectrum. 

The obtained results are illustrated in Fig. 5 and 6. 
The bulk of the gamma-ray emission from this region 
is explained by means of the above described compo- 
nents, but a residual emission is left. The systematic 
uncertainty of the effective area of the LAT is 10% at 
100 MeV, decreasing to 5% at 560 MeV and increas- 
ing to 20% at 10 GeV. This uncertainty propagates 
to the model predictions and should be considered in 
interpreting the residual spectrum in Fig 6. 
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V. CONCLUSIONS 

We have reported a preliminary analysis of the 
Fermi/LAT observations of the Galactic Center. This 
analysis was focused on the indirect search for DM. 
From this preliminary analysis the following conclu- 
sion are obtained: 

• Any attempt to disentangle a potential DM sig- 
nal from the galactic center region requires deep 
understanding of the conventional astrophysics 
background; 

• The bulk of the gamma-ray emission from the 
GC region is explained with the detected sources 
and the Galactic Diffuse emission model, but ... 

• ... a residual gamma-ray emission is left, not 
accounted for by the above models. 

Improved modelling of the Galactic diffuse model 
as well as the potential contribution from other as- 
trophysical sources (for instance unresolved point 
sources) could provide a better description of the data. 
Analyses are under-way to investigate these possibili- 
ties. 
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